Gravitational collapse of Type II fluid in higher dimensional space-times 
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We find the general solution of the Einstein equation for spherically symmetric collapse of Type 
II fluid (null strange quark fluid) in higher dimensions. It turns out that the nakedness and curva- 
ture strength of the shell focusing singularities carry over to higher dimensions. However, there is 
shrinkage of the initial data space for a naked singularity of the Vaidya collapse due to the presence 
of strange quark matter. 
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I. INTRODUCTION 

The gravitational collapse of sufficiently massive star, 
under fairly general conditions, will result into a singu- 
larity is the fact established by the singularity theorems 
of Hawking and Ellis |Q. However these theorems do 
not indicate whether the singularity will be hidden be- 
hind the event horizon or whether it will be visible to an 
outside observer. The singularity will not be visible is 
indicated in the relativistic literature by the phrase cos- 
mic censorship conjecture (CCC). It will not be visible to 
an observer outside the event horizon is the weak form 
while it will not be visible to any observer, even the one 
who is sitting on the collapsing star, is the strong form of 
the conjecture. This conjecture was articulated by Pen- 
rose § (see [| for reviews). The weak form essentially 
states that gravitational collapse from a regular initial 
data never creates the space-time singularity visible to 
distant observers, i.e., any singularity that forms must 
be hidden within a black hole. For the strong form on 
the other hand, nothing is supposed to emanate out of 
the singularity and hence it is not visible to anyone even 
if one is infinitesimally close to it. 

There is no general theory of nature and visibility of 
singularities. There do exist a number of exact solutions 
of the Einstein equation which admit depending upon 
the initial data black holes (BH) or naked singularities 
(NS) (see, e.g., Joshi ||). In particular the Vaidya so- 
lution Q is extensively used to show that the end state 
of collapse for a regular initial data results into a naked 
singularity. In the absence of a general result, the study 
of various examples of collapse becomes pertinent to ex- 
amine the validity of the phenomenon of occurrence of 
naked singularities. The two most investigated cases are 
that of inhomogeneous dust and the Vaidya null fluid. A 
more general matter distribution should be considered. 



Undoubtedly, as singularity is approached, matter would 
be in highly dense and exotic state. The strange quark 
matter is the densest state of matter known which is pro- 
duced in quark-hadron phase transition in the early Uni- 
verse or at ultra high energy neutron stars converting 
into strange quark matter stars gj. It would therefore 
be interesting to study the collapse of the mixture of the 
Vaidya fluid with the strange matter (SQM), which is 
a Type II fluid pL This has been studied in the usual 
4-dimensional space-time and here we wish to study 
it in higher dimensions. The presence of SQM tends to 
shrink the initial data window leading to NS. As is ex- 
pected, the increase in dimensions would also tend to en- 
hance the shrinkage because of gravity getting stronger 
in higher dimensions H. 

Recent developments in string theory indicate that 
gravity may be truly higher-dimcnsional(HD) interac- 
tion, becoming effectively AD at lower energies. The 
main question is how to bring this theoretical develop- 
ment on the anvil of observations? One of the ways is to 
study its effects in astrophysical settings. That is where 
the phenomenon of gravitational collapse in HD attains 
relevance. In this study we would be finding the effects 
of HD as well as SQM on the end product of collapse 
of null SQM. For this, we first obtain the exact solution 
describing collapse of null SQM fluid sphere, and then 
investigate their effect in the context of CCC. 

In section II, we obtain the generalized Vaidya solu- 
tion for null quark fluid in higher dimensional spherically 
symmetric space-time. For definiteness we shall refer this 
metric as HD generalized Vaidya metric. We adapt this 
solution to study end state of collapse in section III . Fi- 
nally we shall conclude with a discussion on the effect of 
SQM on HD Vaidya collapse. 



II. TYPE II FLUID COLLAPSE IN HIGHER 
DIMENSIONS 
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Let us begin with the general (n + 2)-dimensional 
spherically symmetric space-time, in advanced Edding- 
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ton time coordinate v, described by the metric [|[ 9: 



ds 2 
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where < r < oo is the proper radial coordinate, — oo < 
v ^ oo is an avanced time coordiante, and 
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is the metric on n-sphere in spherical polar coordinates. 
It proves useful to introduce a local mass function m(v, r) 
defined by / = 1 — 2m jr [|l0| . Here m(v,r) is an arbi- 
trary function of advanced time v and radial coordinate 
r. When m = m(t>), it is the Vaidya solution in higher 
dimensions j?], g|. The usual Vaidya solution in 4-space- 
time follows for m = m(v) and n — 2. 

We wish to find the general solution of the Einstein 
equation for the matter field given by Eq. (||) for the met- 
ric (Q), which contains two arbitrary functions. Hence it 
is general in the retarded coordinates used. It is the field 
equation G\ — that leads as in the 4D case as well as 
in the present HD case to e^"' r ) = g(v) ||. This could 
be absorbed by writing dv = g(v)dv. Hence, without loss 
of generality, the metric (2) takes the form, 



ds 2 



1 - 



2m(v, r) 



(n - l)r n ~ l 



dv 2 + 2dvdr + r 2 dfll (3) 



The energy-momentum tensor for the Type II (null 
SQM) fluid is given by 0, |§, 

Tab = plJb + (P+P) {laVb + IbVa) + P9ab (4) 



l a = 6°, n a 
1 = 51, n a = 



2m(v, r) 
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l a l a = n a n a - l a n a 



-1. 



(5) 



Here p, p are the strange quark matter energy density 
and thermodynamic pressure while /i is the energy den- 
sity of the Vaidya null radiation. The null vector l a is 
a double null eigenvector of T a b- Physically occurring 
distribution is null radiation flowing in the radial direc- 
tion corresponding to p = p = 0, the Vaidya space-time 
of radiating star. When p = 0, T ab reduces to degen- 
erate Type I fluid B, and it represents string dust for 
p = = p. The energy condition for such a distribution 
are as follows 0: 

(i) the weak and strong energy conditions, 



p>0, p>0, p>0 



(6) 



(ii) the dominant energy condition, 
p>0p>p>0. 



(7) 



By choosing properly the mass function m(v,r), it is pos- 
sible to satisfy the energy conditions. In the particular 
case of m = m(v), the energy condition reduces to p > 0. 
The field equations are S, 



8irp = 
8np = 



(n — l)r n 
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(8b) 
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Here and in what follows dash and dot denote derivative 
respectively d/dr and d/dv. The part pl a h of T a b in (0) 
is the component of matter field that moves along the null 
hypersurface v = const.. In particular when p = p = 
we have the Vaidya solution in higher dimensions. Thus 
the distribution in (Q) represents Vaidya radiating star 
in Type II fluid universe in higher dimensions. Note that 
when p = 0, it is static perfect fluid space-time, which 
will yield regular space-time only when m oc r n+1 which 
is the de Sitter space. As in the 4.D case, the equation 
of state for the zero temperature quark matter is taken 
& 
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where k > is a constant. Inserting the above equation 
in the field equations, we find that 



32nB 
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If we choose the function m(v,r) such that 
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where mo(v, r) is an unknown function, then we find that 
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n(k + 1) ' 
It admits the general solution 

too = S{v)r- n/k+1 +M(v). (13) 
Now, one can easily calculate m, p, p and p, which are 

A 
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TABLE I: Variation of critical parameter a c and yo with D in HD Vaidya and HD null SQM Collapse 



D = n- 



HD Vaidya Collapse 



HD null SQM Collapse (J3 = 0.001, k = 3) 



Double Roots(yo) 



Double Roots(yo) 



0.125 
0.037037 
0.0131836 
0.00512 



4 
3 

2.66667 
2.5 



0.12437012771 
0.036037037037 

0.01179754129 
0.0032850449803 



4.00339 

3.0 
2.65859 
2.47188 



This is the general solution of collapsing null strange 
quark fluid in higher dimensions. The metric explicitly 
reads as 



ds 1 



2M{v) 



2S(v) 



A 



n 2 — 1 



(n - l)^™- 1 ) (n - l) r [(«-2)fc+n]/fc 
dv 2 + 2dvdr + r 2 dfl 2 . (18) 



Clearly, all the energy conditions would be satisfied for 
n > 2 because it would ensure p > and p > 0, while 
fi > would be taken care of when we choose the mass 
functions for both null radiation and SQM. The case, 
S(v) = 0, corresponds to the HD Vaidya de Sitter solu- 
tion 0, D, whereas the Vaidya solution is obtained by 
taking S(v) = B = 0. The solutions in Refs. [|[ |l3| can 
be recovered by setting n — 2. In the next section, we 
shall adapt the solution obtained to study the end state 
of the collapse for k = 3 in Eq. (^) , for which the matter 
is quark plasma. 



III. EXISTENCE OF A NAKED SINGULARITY 



Clearly the above differential equation has a singular- 
ity at r = 0, v — 0. We note that a singularity forms 
when a shell hits the origin at r — 0, where the den- 
sity diverges. If the singularity is at least locally naked 
(for brevity we have addressed it as naked in this pa- 
per), there must be a light ray coming out from it. The 
critical ray is the Cauchy horizon defined as the first 
outgoing radial null geodesic emerging from the singu- 
larity. Therefore, by investigating the behavior of radial 
null geodesies near the singularity, one can find out if 
outgoing null geodesies meet the singularity in the past. 
In order to determine the nature of the limiting value 
of y at r = 0, v — on a singular geodesic, we let 
yo = lim r _ >0 „->o V = lim T -, 0v ^ v/r. Using Eq. (|2C 
and L'HopitaPs rule we get 



yo = lim y = lim — = lim — — = 



r^0 v^0 



lim 

r— >0 v^Q 1 — Oty Tl 



i — >6 v— »o r r—>6 v—>o dr 
2 



P y n(k+l)/k 



A r 2 



(21) 



The physical scenario here is that of a radial influx 
of null fluid in an initially empty region of the higher 
dimensional non-flat but empty space. The first shell 
arrives at r — at time v = and the final at v = T. A 
central singularity of growing mass is developed at r = 0. 
For v < we have m(v) — 0, i.e., the higher dimensional 
de Sitter like space, and for v > T, m(v) — 0, m(v) 
is positive definite. The metric for v = to v = T is 
the higher dimensional generalized Vaidya metric, and 
for v > T we have the higher dimensional Schwarzschild 
solution. For the HD Vaidya region, we choose, 



2M(v) = a(n- l)u n_1 (a > 0), (19a) 
2S(v) = (3(n - l) w K»- 2 )*+»]/* ((3 > 0), (19b) 

Radial null geodesies of the metric (|J) must satisfy 

2M(v) 2S(v) A 



dr 
dv 



1 - 



(n — l)r( n_1 ) (n - i) r [(™-2)fc+«]/fc 



(20) 



which can be rearranged as 



[(n-2)k+ra]/fc+l 



+ ayZ - y Q + 2 = 



(22) 



This algebraic equation is the key equation which gov- 
erns the behavior of the tangent vector near the singular 
point. 



If the singularity is naked, Eq. fl22| ) must have one or 
more positive roots yo, i.e., at least one outgoing geodesic 
that will terminate in the past at the singularity. Hence 
in the absence of positive roots, the collapse will always 
—? ,kad to a black hole. Thus, the occurrence of positive 
1 repots implies that the strong CCC is violated, though 
not necessarily the weak one. 



The limit (3 — > corresponds to HD Vaidya collapse 0, and in that case, Eq. ( p2|) admits positive roots for a < 
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TABLE II: Variation of (3 T with D. For f3 > (3 T , the end state 
of collapse is a black hole for all a (a^ = 10~ 6 ) and k — 3 



D = n + 2 


Critical Value (3t 


Equal Roots yo 


4 


0.20519540362596 


4.99999 


5 


0.037036037037 


3.0 


6 


0.009547734436 


2.6 


7 


0.002807259781 


2.42857 


8 


0.00088477111 


2.33334 



a v c . Hence singularities are naked for a € (0,a c ], and 
black holes form otherwise. Thus a — a v c is the critical 
value at which the transition occurs, the end state of 
collapse switches from naked singularities to black holes. 
This is in agreement with earlier work (see, e.g., 

We want to investigate the changes, at least qualita- 
tively, in the above picture of HD Vaidya collapse in the 
presence of a SQM. It is verified numerically that as n in- 
creases, the threshold value [3t decreases (see Table ||). 
For (3 > (3t, it is BH for all a. Further, the critical 
parameter a s c for HD SQM is smaller than correspond- 
ing critical parameter a v c for HD Vaidya , i.e, the naked 
singularity spectrum of the collapsing HD Vaidya region 
gets covered and hence the initial data space (0, a c ] for 
a naked singularity in HD Vaidya collapse contracts due 
to the presence of SQM. Thus, at least at a qualitative 
level, the presence of SQM facilitates the formation of 
black holes. However, it is interesting to see that for 
each (3, there exists a a s c such that singularities are al- 
ways naked for all a G (0, a s c ]. That is for each (3 there 
exists a non zero measure set of a values giving rise to 
NS and consequently violating CCC (see Tables Q, 



Other Cases 



a. Case k = n. 
This corresponds to 



M = 



87t(n - l)r" 



M(v) + S(v) 



iA 



-P = 



8n(n - 1) 
(23) 



and the algebraic equation takes the form 
(a + 0)y$ - y + 2 = 



(24) 



The metric in this case takes the form of the HD Vaidya 
- de Sitter metric, i.e., the HD collapse of null fluid in an 
expanding de Sitter back ground where A is generated by 
the bag constant. The singularity is naked for 



(a + (3) < A c = 



1 



1 



2n 



The equal roots at Ac is y — 2n/(n— 1). It is remarkable 
to note that both critical parameters and tangents to 



outgoing geodesies are dependent on dimension of space- 
time. Further, we note that j/o is bounded below by value 
2, yo — > 2 as A — > or D — » oo. 

b. Case k — > oo. 
We have 



167r(n — l)r 2 

n0 _ y n-2 



[ay 



n-2 



n— 3i 



p = 0, 



(25) 



I6irr 2 ' 

and the algebraic equation becomes: 

ayo+Pvr 1 -2/o + 2 = 0, 

For n = 2, this would admit a positive root for a < 
l/8(/3— l) 2 , giving the range for naked singularity as ob- 
tained in [ |14[ . This is simply the null fluid collapse in 
the background of constant potential which is character- 
ized by Tg = Tl — const. /r 2 , as is the case in Eq.([25|) 
above. The condition for naked singularity modifies to 
/3 2 < 4a + 108a 2 + 36a^ + (3 2 when n = 3. 
c. Other k 

We also note that as k increases, so does the threshold 
value 0t- However, to conserve the space we shall not 
present the details. 

IV. DISCUSSION 

Motivated by the development of superstring and other 
field theories, there is considerable interest in the models 
with extra dimensions from viewpoint of both cosmology 
[ p| and gravitational collapse |l6[ ]. In this paper, We 
have shown that the 4D spherically symmetric solution 
describing null SQM fluid go over to (n + 2)-dimensional 
spherically symmetric solution and essentially retaining 
its physical behavior. The solution obtained includes 
many previously known solutions to Einstein equations, 
such as, HD Vaidya de Sitter solution. We have also used 
this solution to study the end state of collapsing star and 
showed that there exists a regular initial data which leads 
to naked singularity. Higher dimensions tend to shrink 
the naked singularity initial data space alternatively en- 
largement of black hole initial data space, of correspond- 
ing AD collapse. We have also studied the effect of SQM 
on Vaidya collapse (both in AD and HD). The naked sin- 
gularity spectrum of Vaidya collapse shrinks further with 
introduction of SQM and thus atleast qualitatively SQM 
favors black hole in comparison to naked singularities. 
It is also interesting to see that there exists threshold 
value (3t, which decreases with increase in n, such that 
for (3 > (3t, the end state of the collapse is always black 
hole for all a. However the threshold touches zero only as 
n — » oo. That means for finite dimensions, there would 
always be a window howsoever narrow in the initial data 
leading to NS. However for sufficiently large n, it would 
practically be black hole and CCH would be obeyed. This 
case includes several previous cases of spherical collapse 
in AD and HD. 
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